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Abstract. Results of photometric surveys have brought to light the existence of a population of giant planets 
orbiting their host stars even closer than the hot Jupiters (HJ), with orbital periods below 3 days. The reason 
why radial velocity surveys were not able to detect these very-hot Jupiters (VHJ) is under discussion. A possible 
explanation is that these close-in planets are short-lived, being evaporated on short time-scales due to UV flux of 
their host stars. In this case, stars hosting transiting VHJ planets would be systematically younger than those in 
the radial velocity sample. We have used the UVES spectrograph (VLT-UT2 telescope) to obtain high resolution 
spectra of 5 faint stars hosting transiting planets, namely, OGLE-TR-10, 56, 111, 113 and TrES-1. Previously 
obtained CORALIE spectra of HD189733, and published data on the other transiting planet-hosts were also used. 
The immediate objective is to estimate ages via Li abundances, using the Ca II activity-age relation, and from the 
analysis of the stellar rotational velocity. For the stars for which we have spectra, Li abundances were computed 
as in Israelian et al. (2004) using the stellar parameters derived in Santos et al. (2006). The chromospheric activity 
index Sus was built as the ratio of the flux within the core of the Ca II H & K lines and the flux in two nearby 
continuum regions. The index Sus was calibrated to Mount Wilson index Smw allowing the computation of the 
Ca II H & K corrected for the photospheric contribution. These values were then used to derive the ages by 
means of the Henry et al. (1996) activity-age relation. Bearing in mind the limitations of the ages derived by Li 
abundances, chromospheric activity, and stellar rotational velocities, none of the stars studied in this paper seem 
to be younger than 0.5 Gyr. 

Key words, stars; abundances - planetary systems - planetary systems: formation 

1. Introduction 10 days) it was soon realized that migration has to play 

I 1 1 1 a key role in their formation, since a solid core massive 

Since the discovery of 51Peg b byl Mayor fc Queloz| 119950 ^^^^^-^ ^^^^^^^ ab le to grow at s uch a short 
the number of new extra- solar planets has been rapidly distance from the parent star l|Pollack et alJll99i .With 
growmg. More than about 160 extra-solar planets have increasing number of discoveries, statistical trends are 
been found and yet none of them look like our paradigm ^^-^-^^^ revealing interesting imprints thought to be left 
for planetary formation, i.e., the Solar System. On the migr ation mecha nism which brought these plan- 
contrary, these new discoveries provide evidence that a JUdrv et al.ll20oi . An annoying point when 
large variety of orbital characteristics probably related to ^^^j^^g ^^^^ migration is how to stop it. The pile-up of 
their formation process are actually possible. p^^^^^^ ^^^^ p^^^ ^3-3.5 days in the the orbital period 
Due to the considerable number of giant planets found distribution was understood as an observational clue indi- 
orbiting their host stars at short periods (typically below mating the migration parking point. Although appealing, 

this interpretation was severely challenged by the discov- 

* Based on observations collected at the La SiUa Parana ^^^^^ transiting planets. Some of these new transiting 

Observatory, ESO (Chile) with the CORALIE spectrograph at i , i ,i t ■■ i . .i i 

- 1 r. . . 1 1 TTimr. . irfrnTrTmr, plaucts havc morc thau a Jupiter mass but nevertheless 

the 1.2-m Euler Swiss telescope and UVES at the VLT/UT2 i , • i 

, , • riTc n-\or\ they revolve around their host stars m an even shorter 
Kueyen telescope (Observing run 075.C-0185). 

1 For a continuously updated hst see table at Orbital period than the so-called Hot Jupiters (HJ), typ- 
|http://obswww.unige.ch/Exoplanets, 
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ically Porb ^1 — 2 days. For this reason, they have been 
called Very Hot Jupiters (VHJ). 

Why these planets have not been detected more of- 
ten by the r a dial v elocity surveys is a matter of debate. 
iGaudi et al.l l)2005|) advocate that the non-detection of 
VHJ by the radial velocity surveys is a statistical effect 
due to the much lower frequency of VHJ as compared to 
that of HJ. 

Given their proximity to their host stars, the VHJ may 
experience some degree of evaporation due to the heating 
by stellar UV photons. Vidal-Madjar et al. (2004) showed 
that HD209458b is indeed evaporating, but due to satura- 
tion of the spectral lines, only an upper limit on the mass 
loss of 10^" gs~^ can be derived. According to theoretical 
calculations, the extreme case of an efficient heating of the 
planet atmosphere may lead to an energy-limited evapora- 
tion 1012 gs-i for HD209458b (Lammer et al. 2003), and 
hence to a very efficient evap oration, possibly to the total 
dissociation of some planets ((Lecavelier des Etangs et all 
120041 1 Baraffe et al ]|2004). Conversely, cooling in the at- 
mosphere may well lead to a much more limited mass loss 
close to the lower limit derived by Vidal-Madjar et al. 
llYeiyi2004llHubbard et alJl200,'Tl) . 

Is a catastrophic fate a possible explanation for the 
lack of planets with periods shorter that 3 days? The ex- 
istence of the VHJ is an important constraint to this type 
of scenario. Evaporation scenarios must be able to explain 
the observed evaporation rate of HD 209458b, while at the 
same time accounting for the long-term survival of the 
VHJ on even closer orbits. One possibility is that VHJ 
are not viable in the long term and are fated to catas- 
trophic evaporation. In this case the three observed VHJ 
would simply be too young to have evaporated yet. 

The aim of the present article is to examine the like- 
lihood of the hypothesis that the three VHJ were found 
around parent stars that are very young compared to typi- 
cal planet host stars and to the time scales of evaporation. 
For this purpose, three indicators of young stellar ages 
are c ombined: i) chromospheric activity (e.g. iHenrv et al.l 
19961) ii) Lithium abundances and iii) rotational velocity. 

The manuscript is organized as follows. Our observa- 
tions and data reduction are described in Sec. 2. In Sec. 3 
our index Sus is built and calibrated with respect to the 
index S of Mount Wilson (Smw)- In Sec. 4, we derive the 
Li abundances for the observed stars. Projected rotational 
velocities available in the literature for stars hosting tran- 
siting planets are compiled in Sec. 5. Using the different 
observables we derive ages in Sec. 6 and discuss the im- 
plications of the derived ages for the planetary evolution 
under evaporation in Sec. 7. 

2. Observations 

For 5 of the transiting planet host stars (OGLE-TR-10, 
56, 111, 113 and TrES-1) the observations were carried 
out with the UVES spectrograph at the VLT-UT2 Kueyen 
telescope (program ID 75.C-0185), between April and May 
2005 (in service mode). The Dichroic #1 390-^580 mode 
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Fig. 1. On line and continuum filters used in the compu- 
tation of Sus- The UVES spectrum for TrES-1 (top panel) 
and OGLE-TR-10 are shown in the figure. 

was used. The red arm spectra cover the wavelength do- 
main between 4780 and 6805A, with a gap between 5730 
and 5835A, whereas the data obtained with blue arm cover 
the wavelength domain between 3260 and 4450A. For all 
stars, we adopted slit widths of 0.9 arcsec and 1.1 arcsec, 
which provides a spectral resolution (A/ A A) of the order 
of 50 000 and 40 000 in the red and blue arm, respectively. 

For each exposure, both the blue arm CCD and the 
red mosaic were read in 2x2 bins to reduce the readout 
noise and increase the number of counts in each bin. This 
procedure does not compromise the resolving power, since 
the sampling of the CCD is still higher (by a factor of 2) 
than the instrumental PSF. For the brighter TrES-1, we 
choose a 1x1 binning for the red arm CCD. 

For the OGLE stars, particular attention was given to 
the orientation of the slit given the relative crowdedness 
of the fields. The angle was chosen in each case using the 
images available at the OGLE website^, in order that no 
other star was present in the UVES slit. 

Data reduction was carried out by the ESO staff us- 
ing the UVES pipeline as part of the service mode data 
package. 

^ http://www.astrouw.edu.pl/ ftp/ogle/index. html 
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Finally, for the analysis of HD189733 we have used 
the sam e CORALIE spectrum taken by iBouchv et alJ 
ll2005bt) . 

3. Chromospheric Activity 

Chromospheric ages are commonly derived using the age- 
chromospheric activity relation constructed based on the 
Mont Wilson (MW) c hromospheric flux index (5'Mw)(e.g. 
IVaughan et"al] 119781: iHenrv et al]ll99(i() . In order to use 
MW calibration, one needs first to build a chromospheric 
flux indicator (here called Sjjs) and then calibrate it with 
respect to the MW system (Smw)- This procedure has 
been carried out many times in different papers (e.g . 
ISantos et all 1200(11 IWright et a l."2004': Saffc et al. 2005). 
Here we follow the prescription of .Santos et al.. (,200Q'l . 

The Sus index is defined as {Fh + Fk)/{Ch + Cr) 
where Fjj and Fk are the sum of the fiuxes within a 1-A 
box centered on the Ca H and K lines, respectively. Ch 
and Ck are the total fiux computed as for Fh and Fk but 
for 20-A boxes centered at 4001.067A and 3901.67A , re- 
spectively. In order to calibrate Sjjs as a function of Smw 
a number of common stars in both samples is needed. 
Our original sample observed with UVES contains only 
spectra of stars hosting transiting planets, therefore exter- 
nal data for calibrators is needed. HARPS spectra, kindly 
provided by the HARPS team, were used as calibrators. 
Therefore, as a prior step before computing the flux index, 
the HARPS spectra were smoothed to match the UVES 
blue-arm resolution of 40000 by convolving a gaussian pro- 
file with the appropriate width and the UVES data were 
rebinned to match the HARPS spectral sampling of O.OlA 

Spectra of calibrators and planet host stars were 
brought to the rest frame using the radial velocities given 
by the cross-correlation function fit. The data were then 
trimmed from 3870A up to 4070A and normalized by a 
straight line within this interval. Finally the Sjjs is com- 
puted by adding the relative intensity of the pixels within 
each filter. Figure^shows the spectral region used to com- 
pute Sus- 

One might inquire whether a calibration built with 
data collected with one instrument can be used to cal- 
ibrate a different one. The index S simply measures the 
ratio between the fiux in the core of the Ca II H and K lines 
with respect to a continuum region next to these lines. 
Therefore provided that i) the instrumental responses are 
correctly subtracted within the region containing the Ca 
II lines and the continuum, ii) the fiuxes are computed 
exactly in the same way for both instruments and, iii) 
data have the same sampling and resolution, the indices S 
will refiect solely the stellar contribution regardless of the 
instrument used to collect the data. 

Most of the instrumental contribution (blaze function, 
pixel-to-pixel variations, CCD response, etc.) is removed 
by the division of the extracted spectrum by the flat flcld 
spectrum. The remaining effects are subtracted by nor- 
malizing the reduced spectrum by a pseudo-continuum. 



Table 1. Smw computed from our chromospheric index 
flux Sus using Eq. ^ and the chromospheric fluxes cor- 
rected from the photospheric contribution (< logi?^^ >). 



Dtar 


D- V 


Smw 


< iOg J^HK > 


TrES-1 


0.848 


0.247 


-4.785 


OGLE- TR- 10 


0.606 


0.192 


-4.804 


OGLE-TR-56 


0.589 


0.120 


-5.358 


OGLE-TR-111 


0.933 


0.275 


-4.812 


OGLE-TR-113 


1.020 


0.448 


-4.685 


OGLE-TR-132 








HD149026 








HD189733 


0.898 


0.437 


-4.537 


HD209458 


0.563 


0.154" 


-4.988 



"Taken from lWrieht et alJ ^Q^) 



Since the position of the true continuum is not known, 
the normalization process will certainly introduce errors 
in our computed S and therefore in the final ages. In or- 
der to quantify how changes in the continuum are refiected 
in the computed Sus , we compared the Sus values deter- 
mined using the normalized and non-normalized spectra. 
In the worst cases the variations in the Sus were of 5%. 
This error is much smaller than the variations of Smw 
which are typically of 20%. 

As an additional check of the validity of our Sus — 
Smw calibration, we have compared the Smw derived 
using spectra collected with HARPS with those com- 
puted using spectra f rom the UVES POP spectral library 
l|Bagnulo et"aLll2003|) for 6 stars in common. For 5 stars 
whose Smw is below 0.2, the differences in the Smw are 
below 6%. For the active star HD 22049 {Smw = 0.467 
computed from our Sus using Eq. ^| the relative differ- 
ence is 15% which is largely due to its intrinsic variability. 
As discussed in Sec. 6.1, we have conservatively adopted 
an error on Smw of 20% when estimating our final uncer- 
tainty on the derived ages. 

Using the Sjrs com puted as above and the Smw given 
bv lHenrv et alJ l|l996l) we carried out a least-squares linear 
fit to the data which yields the relation: 

Sus ^ 0.06111 X Smw - 0.00341 (1) 

Figure 121 shows our best- fit. The RMS of the residuals 
of the fit (cTfit) is 8 X 10^*. As pointed out bv^ Santos et alJ 
(200(11 . a fit includes not only uncertainties in the calibra- 
tion procedure but also intrinsic stellar variation. Eq. ^ 
can now be used to compute the Smw- 

The Ca II H and K flux corrected for the photospheric 
flux (R 'u ^ — Rhk —Rphot) is computed usine lNoves et alJ 
l|l984l) prescription. The derived S]\iw yielded by Eq. ^ 
and the corrected chromospheric flux R'^k are given in 
Tabled We note that in this procedure for the OGLE stars 
and TrES-1 we have used B-V colors that were derived by 
in verting the B - V vs. (Teff,[Fe/H]) calibration presented 
in lSantos et"al] l)2004|) . 

No spectra in the Call-line regions were available for 
HD149026 and OGLE-TR-132, and no values for the S 
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Table 3. Estimated V sini for the transiting planet host 
stars studied in this paper. 



Fig. 2. Sus compared to the Smw values for a set of 
calibrators. The solid line is the best least-squared linear 
fit. 

index are available in the literature for these stars either. 
For HD209458 we u sed the value Smw = 0.154 derived by 
IWrierht et al.l ll2 004l). the eff ective temperature and metal- 
licitv of lSantos et all l|2nn4 to compute the R'jjx given in 
Table HI 

4. Li abundances 



Lithiu m abundances were obtained as in llsraelian et al.l 
l)2004|) . from a LTE analysis using a recent ver sion of 
the radiative transfer code MOOG llSnedenlll97i) and a 
grid of lKuruc3 l)l993j) plane-parallel model atmospheres. 
Equivalent widths (EWs) for the Li line near 6707. SA (or 
an upper limit for these) were measured using the IRAF 
splot routine inside the echelle package. This procedure 
gives good results since in our spectra the Li line was never 
strongly blended with any other line in the same spectral 
region. 

The measured EWs (or upper limits for these) were 
then used to derive the Li abundances. The results are 
pr esented in Tabled W e used the stellar parameters listed 
in lSantos et al.l l|2006^ . also listed in Table[21for complete- 
ness. For more details on the derivation of L i abundances 
we refer the reader to llsraelian et alJ ^^2004^ . 

For HD2 09458 we have taken the Li abundance from 
the study of llsraelian et al.l l|2004|) . No Li abundances are 
available for both HD149026 and OGLE-TR-132. 

5. Projected rotational velocity 

For most of our candidates the projected rotational 
velocities and the respective erro rs listed in TableC 



were takeri frorn the literature dOueloz et alJ l200flt 
Pont et all l2004t iBouchv et all I2004r 2005at ISato et all 



ront et ai.ii:^nn4it moucnv et ai.i |zyu4, i^iato et ai 

2005HLaughlin et alJ2005l) . For OGLE-TR-10 and OGLE 



TR-56 the projected rotational velocities were derived 
from HARPS spectra using a calibration based on the 
Cross -Correlation Function width (see e.g. ISantos et alJ 

I2nn2|) . 



Star 


\/rr,f fkm S ^1 


Reference 


TrES-1 


1 ns+0 3 


T.aiicrhlin et al ('2005'! 


OGLE-TR-10 


7.0±1.0 


HARPS CCF 


OGLE-TR-56 


3.2±1.0 


HARPS CCF 


OGLE-TR-111 


<5 


Pont et al. f 20041 


OGLE-TR-113 


<5 


Bouchv et al. f 20041 


OGLE-TR-132 


<5 


Bouchv et al. f 20041 


HD 149026 


6.0±0.5 


Sato et al. (2005) 


HD189733 


3.5±1.0 


BouchY_et_al. f2005al 


HD209458 


3.75±1.25 


Oueloz et al. (2000) 



We note that in our case the V sini values corre- 
spond to the real equatorial rotational velocity, as long 
as the rotational axis of the star is perpendicular to the 
orbital plane of the transiting planet. This hypothesis 
is supported by observations of the Rossiter-McLauglin 
effect on H D209458 and HD189733 ijOueloz et al.ll2000l: 
iBouchv et al...2005a^. 



6. Age of stars hosting Transiting planets 

6.1. Chromospheric ages 

The formulae used in deriving stellar ages from chromo- 
spheri c activity indices are summarized in IWright et all 
l|2004|) . A typical projection of the age-activity relation is 
shown in Fig. O The relation between activity and age 
is steep for young ages and gets flatter for older ages. 
Additionally, stellar activity follows long-term cycles, and 
the use of instantaneous chromospheric flux instead of the 
flux averaged over an entire magnetic cycle leads to an over 
or under estimation of age depending on which moment 
of its cycle the star is observed. The classical example be- 
ing the case of the Sun whose logi?^^ varied from -4.75 
to -5.10 during the "Maunder Minimum", cor responding 
to ag es of 8.0 and 2.2 Gyr, respectively fe.g. iHenrv et alJ 
Il99j). 

IPace fc Pascuinil l)2004j) have studied empirically the 
reliability of chromospheric ages by measuring the chro- 
mospheric activity level in dwarfs belonging to 5 different 
open clusters with ages spanning from 0.6 Gyr up to 4.5 
Gyr. They found that after a strong decrease in chromo- 
spheric activity taking place in solar-type stars between 
the Hyades and the IC 4651 age (i.e., 0.6 Gyr and 1.7 Gyr), 
activity remains virtually constant for more than 3 Gyr. 
Hence an activity-age relation holds up to only about 2.0 
Gyr. Therefore, as cautioned by many authors, chromo- 
spheric activity indices are good indicators of youth, but 
cannot yield reliable age estimates beyond ~ 2 — 3 Gyr - 
as illustrated graphically on Fig|ni 

Consequently, since all stars in our sample show low 
levels of activity, we only derive lower age limits rather 
than speciflc age estimates from chromospheric activity. 
We used the R'jjx values listed in Tabled and the call- 



Melo et al.: On the Age of Stars Harboring Transiting Planets 
Table 2. Stellar parameters and Li abundances for all known transiting planets. See text for references. 
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Ttii« 


JJdJJci 


OGLE-TR-56 


6119±62 


4.21±0.19 


1.48±0.11 


0.25±0.08 


Santos et al. (2006) 


2.7±0.1 


This 


paper 


OGLE-TR-111 


5044±83 


4.51±0.36 


1.14±0.10 


0.19±0.07 


Santos et al. (2006) 


< 0.5 


This 


paper 


OGLE-TR-113 


4804±106 


4.52±0.26 


0.90±0.18 


0.15±0.10 


Santos et al. f2006) 


< 0.2 


This 


paper 


OGLE-TR-132 


6411±179 


4.86±0.14 


1.46±0.36 


0.43±0.18 


Bouchv et al. (2004) 








HD 149026 


6147±50 


4.26±0.07 




0.36±0.05 


Sato et al. f2005~l 








HD189733 


5050±50 


4.53±0.14 


0.95±0.07 


-0.03±0.04 


This paper 


< -0.1 


This 


paper 


HD209458 


6117±26 


4.48±0.08 


1.40±0.06 


0.02±0.03 


Santos et al. (2004) 


2.7±0.1 


Israehan et al. (200A) 



0.6 



'OA 



0.2 
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ISestito fc RandichlEii . Due to its fragility, Li is quickly 
destroyed during the pre-main sequence inside solar-mass 
stars which are (almost) fully convective at the begin- 
ning of their lives. Li depletion also occurs on the main- 
sequence altho ugh its causes are still under strong debate 
(see review bv lDeliv annis et al.ll2000|) . For this reason Li 
is widely used a youth indicator fe.g. lMartinlll997j) . 

For stars hotter than 6000 K (i.e., earlier than G), 
the convective zones are not deep enough to bring Li to 
warmer regi ons thus Li is not depleted ( with exception of 
the Li gap; iBoesgaard fc TriT)iccolll983) . Hence Li abun- 
dance cannot be used as an age indicator. This is the case 
for OGLE-TR-10 and OGLE-TR-56. 



The recent results of lSestito fc R.andichI ^ 



seem to 

indicate that Li-depletion occurs by quick episodes rather 
than at a monotonic pace. This could imply that Li abun- 
dances can only yield lower limits for the ages. 



Fig. 3. The projection of the age-activity relation for a 
B — V = 0.63 star is shown by the circles and solid line. 
The shaded area around the solid line is the 95% confi- 
dence region for the true stellar age computed assuming a 
20% variability level on the Sjjs- Due to the the flatness 
of the activity-age calibration it is clearly seen that af- 
ter 2-3Gyr only a minimum age can be computed by this 
method. 



Table 4. Age constraints yielded by each method. The 
95%-confidence interval (2-sigma) for the chromospheric 
age computed using the mean value of log R'j^j^ is shown 
along with and the age constraint from Li abundances. 
Whenever available age estimates found in the literature 
are also given. 



bration of lHenrv et al.l l(l996l) . We calculated the 2-sigma 
lower limits on the ages by including as far as possible 
the sources of error on all intermediate steps between the 
computation of the Sjjs and the final ages and allowed for 
an intrinsic dispersion of 20% for the activity level of stars 
with identical mass and age. 

6.2. Li ages 

Li ages given in Tabled] were derived by comparison with 
abundance curves of Li as a functio n of effective tem- 
perature for clusters of different ages Ij.Teffries et al...2002; 



Star 


agecaii 


ageti 


age 


Ref 


2 


— a lower lim. 




other 


other 


TrES-1 


> 1.1 


>0.6 


2.5 ± 1 


1 


OGLE-TR-10 


> 1.1 








OGLE-TR-56 


> 2 




3± 1 


2 


OGLE-TR-111 


> 1.1 


>0.6 






OGLE-TR-113 


> 0.7 


>0.6 






OGLE-TR-132 










HD 149026 






2 ±0.8 


3 


HD 189733 


> 0.5 


>0.6 






HD209458 


> 2 




4.5 


4 



l-'Sozzetti et alJ 120041); 2-'Sasselov 
3 -.Sato et aL (.2005.1 : 4 -.Mazch ct alj 
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6.3. Vsini ages 

During the prc-main sequence phase, angular momentum 
evolution is mainly driven by the magnetic cou pling be- 
tween the star and its disk (e.g. iKonigll llQQlj) . Due to 
differences in the disk-locking time-scales, solar type stars 
arrive on the zero age main sequence presenting a large 
spread in their rotation rates. From this point on, any 
angular momentum evolution is dictated by stellar winds 
whose intensity is itself a functi on of magnetic activity 
and rotation (e.g. lKawaleilll988l) . This creates a regula- 
tion mechanism leading to a steady decrease and to the 
eventual convergence of the rotation rate at about IGyr. 

Therefore the use the rotation rate as a qualitative 
age diagnostic is limited. Neverthel ess, comparing^ our \^rot 
given m Table El with Figure 2 of iBouviej l )l997(l we see 
that despite the large dispersion in the points, the F and 
G stars with Vrot < 5 km s"""^ are likely to be older than 
the Hyades (~ 600 Myr), whereas for the cooler stars a 
Vrot of 5 km is still compatible with the age of the 
Pleiades 100 Myr). Similar conclusions can be drawn 
for OGLE-TR-10 or HD149026 whose Vrot = 7 and 6 km 
s""'^, respectively, can still be found (although unlikely) 
among the Pleiades members. 

The utility of the Vrot values to derive stellar ages is 
limited, but our results show that globally speaking our 
sample stars, and in particular those orbited by VHJ, are 
ver y unlikely to be younger than the Pleiades (see Figure 
2 of lRouvieT<IT997|) . 

A summary of the age constraints yielded by each 
method is given in Table 0] As dicussed in Sec. 6.1, any 
attempt to assign a precise age based on the chromo- 
spheric activity is meaningless due to the flatness of the 
age-activity calibration beyond 2-3Gyr. Therefore, only 
the lower limit (2-sigma level) from chromospheric activ- 
ity and the age constraint from Li abundances are given 
rather than a precise age. Age estimates found in the liter- 
ature, quoted in Tabled are consistent with these values. 

7. Constraining planetary evaporation rates from 
stellar ages? 

The lower limits given in Table 4 allow us to state that 
the stars harboring transiting planets studied here are not 
younger than 0.5 Gyr and for most of them not younger 
than 1 Gyr at a 2-sigma level. Is this compatible with a 

very efRcient evaporation of these o bjects? 

Efficient evaporation scenarios ('Lammer et al.' '2003': 
iLecavelier des Etangs et al 2004.: Baraffe et al 2004') im- 
ply that a fraction of the close-in planets do not survive, at 
least not as gaseous giants. Therefore those planets seen in 
present-day transit surveys correspond to the lucky ones 
that were initially massive enough to escape (so far) evap- 
oration. Our age determinations allow us in principle to 
constrain the magnitude of the evaporation: stars harbor- 
ing very close-in planets (and in particular the transiting 
planets discovered so far) should be on average younger if 
an efficient evaporation indeed took place. For very large 



evaporation rates, old stars should have lost their close-in 
planetary companions, except for the rare ones with very 
large initial masses. 

In order to estimate how the mean age of stars 
harboring very close-in planets depends on the 
magnitude of the planetary mass-loss, we proceed 
as follows. First, it is assumed that the initial 
masses of the planets are set according to a fixed 
planetary mass function (see below). Monte-Carlo 
approach is then used to draw a large number 
of star-planet systems for which we calculate the 
fraction of planets that have evaporated away as 
a function of stellar age and evaporation rate. By 
re-normalizing this fraction we derive a distribu- 
tion of the age of stars with planets as a function 
of the assumed planetary evaporation rate. 

The process described above is carried out as- 
suming the following hypotheses: 

1. The initial planetary mass function for very- 
close in planets is the same as for planets 
observed by radial velocimetry with periods 
smaller than 365 days (a value chosen to avoid 
biases for planets with masses > O.SMjup, and 
which is found to affect the results only mod- 
erately). The sample used is taken from J. 
Schneider's web page (www.obspm.fr/planets) 
and contains 99 planets. 

2. Planets disappear from the sample when their mass be- 
comes smaller tha n 0.5M.T„n, either b ecause of a run- 
away evaporation l)Baraffe et al.ll20o'3) . or because the 
presence of a core makes them too small to be detected 
by the transit method fe.g. .Guillot .20051) . 

3. The stellar formation rate (including stars with plan- 
ets) has been constant for the past 10 Gyr. 

With the present understanding of planet formation 
and migration, our first hypothesis should be conservative 
(i.e. allow for higher evaporation rates), as more massive 
planets are expected to migra te less rapid ly and efficiently 
to be come If J or VHJ planets l|Udrv et al.i.2003 : Ida fc Lini 
I2OO4I e.g.). This is also the case of the third hypothesis, as 
the progressive enrichment of the Galaxy in metals by stel- 
lar nucleosynthesis implies that stars with planets should 
be on average younger. More detailed studies are needed 
to address more precisely these points and assess their ef- 
fect on the results. 

Figure^ shows the resulting age distribution of stars 
with planets for different planetary evaporation rates. 
For the most extreme case of a IMjup/Gyr mass loss, 
the mean of the age distribution is lowered to only 2. 5 
Gyr. For comparison, according to ((Baraffe et al.lEoo3) . 
the (maximal) evaporation rates derived from an energy- 
limited approach is almost constant in time and range 
from O.lMjup/Gyr (HD209458b, OGLE-TR-lOb) to al- 
most 0.7Mjup/Gyr (OGLE-TR-132b). 

In addition to the age distribution derived in 
Figure |4l one can now check whether the sample 
of stars with age determinations derived in this 
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Fig. 4. Age distribution of stars with planets, as a func- 
tion of the planetary evaporation rate. We assume that the 
unperturbed sample of stars with planets have a uniform 
age distribution between and 10 Gyr (dashed line). The 
progressively steeper distributions correspond to evapo- 
ration rates of 0.1 (blue), 0.3 (orange) and LOMj^p/Gyr 
(red), respectively. The diamonds represent the mean ages 
of the different distributions. 



paper is consistent with an energy-limited evap- 
oration scenario. Using again a Monte-Carlo ap- 
proach and adopting the evaporation rates derived 
by Baraffe et al. (200 4). the distribution of possi- 
ble ages for each star with transiting planet stud- 
ied is derived. We then compute the number of 
times that these "Monte-Carlo" stellar ages are 
compatible with the inferred minimal ages, both 
in the energy-limited evaporation scenario and in 
the no-evaporation scenario. We found that the 
latter was between 2 and 3 times more likely, i.e. 
that the energy-limited evaporation scenario was, 
in this case, marginally inconsistent with the ob- 
servations. 

Unfortunately, an additional limitation has to 
be considered, namely, the fact that stars with 
young ages (< 0.5 Gyr) are quite active and as 
such are unlikely to be detected as having a planet 
(both by the transit method and by radial ve- 
locimetry). Due to this further constraint and 
given the small sample used in this work, we can- 
not conclude in favor or against the energy-limited 
evaporation scenario. 



8. Conclusions 

We have shown that the analysis of high resolution spec- 
tra of stars with planets could be used to determine lower 
limits on their age using Li abundances, the Ca II activity- 
age relation and from the analysis of the stellar rotational 
velocity. Applied to six stars known to host transiting 
giant planets, OGLE-TR-10, 56, 111, 113, TrES-1, and 



HD189733, we were able to show that none appears to be 
younger than 1 Gyr at the 2-sigma level. 

In principle, this could be used to constrain the 
amount of mass loss experienced by the planets, 
since high evaporation rates should yield smaller 
average ages for the stars with very close-in gi- 
ant planets. Because of the various uncertainties 
and small size of the sample studied here, the re- 
sults remain inconclusive concerning evaporation. 
However, it highlights the importance of a careful 
determination of the ages of stars with planets, 
in particular for the stars harboring very close- 
in planets ( "Pegasids" ) , but also, as a means of 
comparison, for stars with planets on orbits with 
longer-periods. 

We hence recommend to pursue the fine determina- 
tion of the ages of stars with planets, in priority for the 
transiting planets as they ar e progressively bein g detected 
(for a direct application see lGuillot et al.lf2006^ . for non- 
transiting very close-in giant planets (at less than 0.1 AU 
from their star), and then for a sample of stars with plan- 
ets that are orbiting at larger distances. 
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